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Abstract. The primary inversion of the accurately mea- 
sured frequencies of solar oscillations determines the me- 
chanical properties of the Sun, i.e., the sound speed and 
density as a function of solar radius. In order to infer the 
temperature and hydrogen abundance profiles inside the 
Sun, it becomes necessary to use, in addition, the equa- 
tions of thermal equilibrium along with the input physics, 
such as opacities, equation of state and nuclear reaction 
rates. It then becomes possible to estimate the effects of 
uncertainties in opacity and nuclear reaction rates on the 
inferred thermal and composition profiles. The seismically 
determined temperature and hydrogen abundance profiles 
turn out to be close to those in a standard evolutionary so- 
lar model that includes the diffusion of helium and heavy 
elements below the convection zone. The most significant 
departure occurs just below the base of the convection 
zone where the inferred hydrogen abundance profile is 
smoother than that in a standard diffusion model. The re- 
gion just beneath the solar convection zone appears to be 
mixed, a process which could account for the observed low 
lithium abundance in the solar envelope. With a reason- 
able allowance for errors in opacities, the helioseismically 
estimated cross-section for pp nuclear reaction rate turns 
out to be (4.15 ± 0.25) x 1CT 25 MeV barns. 

Key words: Sun: Abundances - Sun: Interior - Sun: Os- 
cillations 



1. Introduction 

The precisely measured frequencies of solar oscillations 
provide us with a unique tool to probe the solar interior 
with sufficient accuracy. These frequencies are primarily 
determined by the dynamical quantities like sound speed, 
density or the adiabatic index of the solar material and 
a primary inversion of the observed frequencies yields the 
soun d speed and density pro files i nside the Sun (Gough 
Gough & Kosovichcv 1990; Gough & Thompson 



et al. 1996; Gough et al. 199E). On the other hand, in order 
to infer the temperature and chemical composition profiles 
additional assumptions regarding the input physics are re- 
quired (S hibaha shi [1993| ; A ntia fc Chitre |1995| ; Shibahashi 
& Takata 1996| ; Kosovichev 1996). Thus, the equations of 
thermal equilibrium enable us to determine the temper- 
ature and hydrogen abundance profiles in the solar inte- 
rior provided the opacities, equation of state and nuclear 
energy generation rates are prescribed. Although the pri- 
mary inversions can yield the sound speed to an accu- 
racy of 0.1%, the opacities and nuclear reaction rates are 
hardly known to comparable accuracy and consequently, 
more systematic errors are introduced in these secondary 
inversions for temperature and chemical composition. 

There are a number of app roach es adopted for sec- 
ondary inversions. Kosovichev (|1996| ) has employed the 
equations of thermal equilibrium to express the changes 
in primary variables (p, Ti) in terms of those in secondary 
variables (Y, Z) and obtained equations connecting the 
frequency differences to variations in abundance profiles. 
It should be noted that modifications in Z profile mainly 
affect the opacities in the solar interior while the equation 
of state and nuclear energy generation rates are affected 
to a much lesser extent. Such a procedure is essentially 
equivalent to finding the Y profile along with the neces- 



sary opacity modifications. Shibahashi and Takata (1996 



hereinafter ST96) adopt the standard opacities and nu- 
clear reaction rates to obtain the temperature and chemi- 
cal abundance profiles using the inverted sound speed pro- 
file. 



1985 



1991; Dziembowski et al. 1994; Antia & Basu 1994a; Basu 
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Antia and Chitre (1995, 1996) set out to estimate the 
central temperature of the Sun. They adopted the inverted 
sound speed and density profiles to obtain the tempera- 
ture (T) and helium abundance (Y) profiles in the solar 
core, but the main difference was that opacity and nu- 
clear reaction rates were not directly employed for this 
purpose. Instead, the T and Y profiles were obtained by 
minimizing the variation in opacities from the standard 
values. The main reason for allowing variations in theo- 
retically determined quantity like opacity rather than the 
'experimentally' inferred seismic sound speed and density 
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was that the uncertainties in opacities are probably larger 
than those in primary inversions. Another advantage of 
this approach was that it enabled us to study the effect of 
changes in opacity and nuclear reaction rates on the in- 
ferred thermal profiles in a straightforward manner. The 
main thrust of the foregoing study was to estimate the 
central temperature of the Sun and the neutrino fluxes. 
On the other hand, Roxburgh (1996) and Antia & Chitre 



( |1997| ) studied various possible abundance profiles with- 
out any bounds on the opacity variations to study the im- 
plication of helioseismic constraints on the solar neutrino 
problem. 

In the present study we extend the earlier work of An- 
tia and Chitre ( |l995| ) to determine the temperature and 
hydrogen abundance profiles throughout the radiative in- 
terior of the Sun and investigate possible uncertainties 
that might exist in the basic nuclear energy generation 
rates and opacities. Recently there has been a claim that 
the pp nuclear reaction rate should be revised upwards by 



a factor of 2.9 (Ivanov et al. 1997) and it would there- 
fore be interesting to test this suggestion helioseismically 
(Degl'Innocenti et al. 1997). Further, in the earlier study 
we had restricted the composition profiles to smooth func- 
tions represented by a low degree polynomial, which con- 
strained the class of admissible solutions. This restriction 
has been relaxed in the present study by using cubic spline 
basis functions to represent the composition profiles. 

The rest of the paper is organized as follows: the inver- 
sion technique employed to obtain the temperature and 
chemical composition profiles is described in Section 2, 
and the results arc set out in Section 3. Our attempts to 
constrain the cross-section for the pp reaction using the 
inverted profiles are outlined in Section 4, while Section 5 
summarizes the conclusions from our study. 

2. The inversion technique 

The sound speed and density profiles inside the Sun are 
inferred from the observed frequencies using a Regularized 
Least Squares technique (Antia 199£). The primary inver- 
sions based on the equations of hydrostatic equilibrium 
along with the adiabatic oscillation equations, however, 
provide only the mechanical variables like pressure, den- 
sity and sound speed. This provides us with the ratio T/p, 
where p is the mean molecular weight. In order to deter- 
mine separately T and p, it becomes necessary to use the 
equations of thermal equilibrium, i.e. 



6Anr 2 aT 3 dT 
3np dr ' 



dL r 
dr 



Anr pe, 



(1) 
(2) 



where L r is the total energy generated within a sphere 
of radius r, a is the Stefan-Boltzmann constant, k is the 
Rosseland mean opacity, p is the density and e is the nu- 
clear energy generation rate per unit mass. In addition, the 



equation of state needs to be adopted to relate the sound 
speed to chemical composition and temperature. Equa- 
tion (Q) is applicable when there is no convective transport 
of energy. This is generally true in the region below the 
outer convection zone, and we have verified that all the 
models considered in this work are stable to convection. 

Since we have only three equations, namely, equations 
(|l]), (j|) and the equation of state to determine the vari- 
ables T, L r and the chemical abundances, it becomes pos- 
sible to determine only one parameter specifying the com- 
position, e.g., the mean molecular weight. Clearly the so- 
lution cannot be unique and therefore, in this work we 
assume the heavy element abundance, Z to be prescribed 
and attempt to determine the hydrogen abundance (X) as 
also the temperature. It should be stressed that Z mainly 
affects the opacity in the solar interior, since the bulk of 
the energy generation takes place through the pp chain. 
Thus, an increase in Z by 20% (which is a reasonable esti- 
mate for possible errors in Z) raises the opacity by about 
8-15%, while the sound speed changes by no more than 
0.2%, and the integrated luminosity changes by at most 
0.4%. It is evident that the dominant effect of a change in 
the Z profile is on the opacities and in this work we do not 
make any attempt to separate the intrinsic errors in opac- 
ity tables from those arising due to uncertainties in the Z 
profiles. The opacity changes could be because of intrinsic 
errors in opacity tables or due to incorrect Z profiles and 
it becomes difficult to separate the two effects. One reason 
to keep the Z profile fixed is that the value of Z/ X in the 
convection zone is known (Grevesse & Noels 1993 ) and the 
change in the interior due to diffusion is not expected to 
be very large, being of the order of 10% (Proffitt 1994| ) or 
eve n less depending on treatment of diffusion (Richard et 
al. 1996 ). We can thus assume the Z profile to be known 
to an accuracy of better than 20%. 

In order to calculate the X profile we express it in 
terms of suitable basis functions, e.g., B-splines, by writing 



i=l 



(3) 



where <fii(r) are the cubic B-spline basis functions based 
on uniformly spaced knots. We use knots with a spacing of 
0.02i? Q , which is found to be adequate to represent the X 
profile to the level of accuracy expected from helioseismic 
inversions. We have tried experiments by increasing the 
number of knots to find that it does not have any signif- 
icant effect on the solution. For a given set of coefficients 
a.i it is possible to calculate X and then the equation of 
state together with the known sound speed and density 
profiles determine the temperature profile. Once the tem- 
perature, density and composition profiles are known we 
can integrate equation (^|) to calculate the luminosity, L r 



Awr 2 pe dr. 



(4) 
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Here the upper limit of integration can be suitably trun- 
cated since there is no significant energy generation in the 
outer layers. It may be noted that this equation is inte- 
grated from the outer boundary in order to prevent the 
errors in primary inversion near the center from contami- 
nating the results in the outer region. With this approach, 
the errors in primary inversion near the center will only af- 
fect the secondary inversion in the central region. We use 
t he nu clear reaction rates from Bahcall & Pinsonneault 
( |l995| , hereinafter BP95) to calculate the nuclear energy 
generation rate e. However, the cross-section for the pp re- 
action which has a dominant influence on nuclear energy 
generation rate in the Sun, has decreased by about 4.5% in 
the last few years (Bahcall |I989 BP95). We therefore, use 
both these values to estimate the influence of uncertain- 
ties in the nuclear reaction rates. Most of the calculations 
have been performed using the older reaction rate for pp 
reaction, as that is found to give the computed luminosity 
closer to the observed value, L©. 

Once the luminosity has been determined in the man- 
ner outlined above, we can substitute all the quantities in 
equation ([[]) to obtain 



R K , = 



2>npL r 



647rrVT 3 ^' 

dr 



(5) 



We use the recent OPAL opacities (Iglesias and Rogers 



1996) to estimate the ratio If the equations of ther- 



mal equilibrium are exactly satisfied, this quantity would 
be unity everywhere, but in general that is not the case. 
The departures of R K from unity is a measure of the extent 
by which the opacity needs to be modified to satisfy the 
equations of thermal equilibrium for a composition profile 
prescribed by equation (|3|), the coefficients in which can 
then be determined by minimizing the required opacity 
modifications. This outlines our prescription for calculat- 
ing the composition and temperature profiles which may 
be regarded as defining a seismic solar model. It may be 
noted that generally with such a procedure the integrated 
luminosity will not turn out to be equal to the observed 
solar luminosity. We can adjust the nuclear energy gener- 
ation rate to obtain the correct luminosity. 

In order to implement this procedure we need to mea- 
sure the required deviation of opacity. It would be sim- 
plest to use a least squares approach, where the integral 
of squared difference is minimized. Thus we minimize the 
quantity 



F = £(ik(r,)-i) a + a( 



L-Lr. 



© n2 



(6) 



where r% are a set of suitable mesh points spanning the 
radiative interior of the Sun, L is the computed luminos- 
ity in the seismic model and a > 0. We generally use a 
mesh with uniform spacing of 0.005i?©, which gives ap- 
proximately 145 points. Thus we determine the X profile 
by choosing the coefficients in equation (^) to minimize 



F. Depending on the value of a, this procedure may also 
be able to yield the correct solar luminosity without ad- 
justing the nuclear reaction rates. Thus, only the opacity 
deviations will be minimized for a = 0, while for large val- 
ues of a, the integrated luminosity can also tend to the so- 
lar luminosity at the expense of larger opacity variations. 
However, in general it is found that an adjustment of the 
luminosity is somewhat difficult because in the process 
the resulting composition profile as well as the required 
opacity modifications may become unacceptable. 

In actual practice the function is linearized about some 
initial guess for the X profile and least squares solution 
is calculated iteratively. Apart from this simple technique 
we have also tried the technique of simulated annealing 
(Vanderbilt & Louie 1984; Press et al. 1993) to obtain the 
nonlinear least squares fit. Since the convergence of sim- 
ulated annealing technique is very slow, after some stage 
we switch to the linearized version to arrive at the final 
solution. 

It is not obvious that such a choice will indeed produce 
the correct X profile since the actual errors in opacity 
may be larger than the minimum estimate, or the opacity 
modification may not be correctly estimated because of 
errors in the primary inversion or in the nuclear energy 
generation rates or in the adopted Z profile. The errors 
in the X profile due to those in the primary inversions 
can be estimated by perturbing the inverted profiles of 
sound speed and density. For this purpose we repeat the 
primary inversions using a perturbed set of frequencies, 
where randomly distributed perturbations with variance 
equal to the quoted errors in frequencies are added to the 
input frequencies. However, it turns out that these errors 
are fairly small as compared to those introduced by un- 
certainties in opacities. For estimating the errors arising 
from uncertainties in opacities we try to determine the X 
profile for various prescribed Z profiles. We use for this 
purpose one of the following three basic Z profiles: (1) a 
homogeneous profile (denoted by HOM) without any dif- 
fusion of heav y elements, (2) a profile including diffusion 
(Proffitt |1994| ) indicated as PROF, and (3) another pro- 
file using a different treatment of diffusion including some 
turbulent mixing just below the base of the convection 
zone (Richard et al. 1996 ) identified as RICH. We scale 
all these profiles to give a prescribed value of Z at the 
solar surface, Z sur { and the inversions are performed for a 
very large range of values for Z SUI { . 

Inside the convection zone we can determine the he- 



lium abundance independently (Gough 1984; Dappcn et 
al. 1988a: Dzicmbowski et al. 1991; Kosovichevet al. 1992; 



Antia & Basu 1994b; Basu & Antia 1995) and the tem- 



perature can then be determined using the inverted sound 
speed. However, in this work we have restricted ourselves 
to the radiative interior only. The estimated value of X 
at the base of the convection zone can be compared with 
the independently estimated values inside the convection 
zone. Once the T and X profiles inside the Sun are known 
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Fig . 1. R elative difference in sound speed and density between the Sun and the Model S of Christensen-Dalsgaard et 
al. ( 1996 ) as inferred using various sets of observed frequencies. The continuous, dotted, short-dashed and long-dashed 
lines respectively, represents the results using GONG data for the months 4-10, month 9, months 4-7 and months 
4-14. While the dot-dashed line represents the results using BBSO data combined with low degree frequencies from 
BiSON. 



it is straightforward to estimate the neutrino fluxes for the 
various solar neutrino experiments. 

The representation (6) takes care of both the modifi- 
cations in the opacity and nuclear energy generation rate. 
Thus, in the limit of a = we essentially recover the 
approach of ST96, for which the required opacity modifi- 
cation is negligible. When we try the least squares solution 
with a = 0, the resulting opacity modifications turn out 
to be generally very small and merely represent numeri- 
cal errors. We have, in fact, verified that the two solutions 
are almost identical. Most of the results in this paper have 
been obtained using this prescription. In our formulation 
by increasing the value of a we can also find solutions 
which satisfy the luminosity constraint at the expense of 
allowing for some opacity variations. This approach would 
be similar to that of Kosovichev (1996) who has estimated 
Z variations, which is essentially equivalent to finding the 
opacity modifications. We thus have the choice of either 
modifying the nuclear energy generation rate to match the 
solar luminosity (a = 0) or to modify the opacities keep- 
ing the nuclear energy generation rate fixed (a » 1) to 
achieve the same purpose. Naturally, by using interme- 
diate values of a we can obtain solutions which require 
modifications of both opacities and nuclear energy gener- 
ation rates by varying amount. Clearly the resulting solu- 
tion is not unique, but all possible solutions may not be 
acceptable since some of them may require unacceptably 
large modifications in opacity or nuclear energy generation 
rates. Of course, if we can get the observed solar luminos- 



ity for a solution with a = 0, then effectively no modifica- 
tion would be required in input microphysics. In fact, for 
some choices of nuclear reaction rates and opacities we do 
find such solutions where no significant modifications in 
microphysics are required. 

It should be stressed that our technique for inferring 
the T and X profiles in the solar interior is absolute, in 
the sense that no reference model is required and the ac- 
tual profiles of T and X are determined directly from the 
sound speed and density profiles. Of course, the density 
and sound speed are determined by using a differential 
technique where the differences are linearized about a ref- 
erence model. 

3. Inverted T and X profiles 

We use data sets of p-mode frequencies from the GONG 
data (Hill et al. 1996| ) and from the Big Bear Solar Obser- 
vatory (BBSO) data (Libbrecht et al. |1990[ ) along wi th the 
low degree modes from BiSON data (Elsworth et al. 1994) 
to infer the sound speed and density profiles using the 
Regularized Least Squares technique (Antia 1996[ ). The re- 
sults obtained using different input frequencies are shown 
in Fig. [I] which displays the relative difference between 
inverted profiles and those in the model S of Christensen- 



Dalsgaard et al. (1996). It can be seen that the results ob- 



tained using different sets of frequencies agree with each 
other to within the estimated errors. The most significant 
difference in the sound speed between the model and the 
Sun, occurs just below the base of the convection zone and 
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Fig. 2. Rela tive difference in sound speed and density between the Sun and the Model S of Christensen-Dalsgaard 
et al. ( 1996 ) as inferred using two different estimates for solar radius. The continuous and dashed lines represent the 
results obtained with estimated radius of 695.99 and 695.78 Mm respectively. The dotted lines show the la error 
limits. 



is very likely on account of the X profile in the Sun be- 
ing smoother than that in the model (Basu & Antia 1994 



Gough et al. 1996; Basu 1997). Apart from this, another 
noteworthy smaller hump occurs around r — 0.2i? Q which 
is opposite in sign to that below the convection zone. It 
is likely that in this region the X profile in the Sun is 
steeper than that in the model. The third minor hump 
around r — 0.05-R© is not very significant and occurs in 
a region where the primary inversions are not likely to be 
in any case very reliable. Similarly the most significant 
difference in the density profile occurs inside the convec- 
tion zone and is probably due to small errors in opacity 
equation of state, surface abundances and/or the depth of 
the convection zone (Basu & Antia 1 9 9 7| ) . 

Recently, it has bee n sugg ested that the standard value 
of solar radius (Allen 1973) needs to be reduced (Antia 



1997 ; Schou et al. 1997 ) and it would be interesting to es- 
timate the effect of error in the solar radius on helioseismic 
inversions. We have shown in Fig. [I] the results obtained 
using the standard value of 695.99 Mm for the solar radius. 
We have also performed inversions with a reduced radius 
of 695.78 Mm and Fig. || compares the results obtained 
with two different values for solar radius using the same 
set of observed frequencies from GONG months 4-10 data. 
It is clear that the small error in solar radius affects the 
inversion results to an extent which is much larger than 
the estimated errors due to those in frequencies. 

Applying the procedure outlined in Section 2 to the 
inverted profiles for sound speed and density shown in 
Fig. [|, we obtain the T and X profiles and the results are 



shown in Fig. ||. All these results have been obtained with 
a = in equation (|^) and using the pp reaction rate from 
Bahcall ( |l989| ). A Z profile including diffusion (RICH) of 
heavy elements with surface value of Z = 0.018 was as- 
sumed for these inversions. Once again it is clear that the 
results obtained using different input frequencies are close 
to one another. The computed luminosity in these seismic 
models turns out to be between 0. 965-0. 993L Q , which is 
roughly consistent with the actual solar luminosity. The 
errors in secondary inversions arising from estimated un- 
certainties in the input frequencies can be calculated with 
a Monte-Carlo simulation. For this purpose we generate 
20 sets of artificial frequency data where randomly dis- 
tributed errors with standard deviation equal to the esti- 
mated errors in observed frequencies are added to every 
input frequency before the primary inversion. The inverted 
sound speed and density profiles are then used to obtain 
the T and X profiles. The standard deviation in these pro- 
files, at a fixed radius, will give an estimate of errors in 
secondary inversions arising from those in input frequen- 
cies. It turns out that the relative errors in the inferred 
values of T are much smaller than those in X, and clearly, 
the temperature gets determined much more reliably than 
the chemical composition in this procedure. 

We can also estimate the influence of error in adopted 
solar radius on the secondary inversions by using sound 
speed and density profiles obtained with different values 
of the radius and these results are shown in Fig. [|. These 
errors can be seen to be comparable to those due to un- 
certainties in frequencies. The properties of these seismic 
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Fig. 3. Relative difference in temperature and absolute difference in the hydrogen abundance X, between the Sun and 
the Model S of Christensen-Dalsgaard et al. ( 1996 ) as inferred using different input frequencies. The line styles are 
same as those in Fig. [|. 
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Fig. 4. Relative difference in temperature a nd abs olute difference in the hydrogen abundance X , between the Sun and 
the Model S of Christensen-Dalsgaard et al. ( 1996 ) as inferred using different estimates for solar radius. The continuous 
and dashed lines represent the results obtained with estimated radius of 695.99 and 695.78 Mm respectively. The dotted 
lines show the la error limits on the results with estimated radius of 695.99 Mm. 



models are summarized in Table |l|, which also gives the es- 
timated errors in each quantity due to those in the GONG 
months 4-10 data. In this table T c is the central temper- 
ature, </>( 37 Cl) the neutrino flux in the Chlorine detector, 
0( 71 Ga) the neutrino flux in the Gallium detector, while 
<fi( 8 B) is the flux of 8 B neutrinos. It can be seen that a 



reduction in radius by 210 km increases the computed lu- 
minosity by O.OO4L0 and the neutrino fluxes are also cor- 
respondingly enhanced by similar amounts. The last line 
in the table gives the results for a static solar model des- 
ignated as INV, which is constructed using the inverted 
X profile. This model will be discussed later in this sec- 
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tion. In the following discussion all the results have been 
obtained using the GONG months 4-10 data with an es- 



timated radius of 695.78 Mm (Antia |1997[ ). 

Apart from this there could be some errors due to un- 
certainties in the equation of state resulting from the use 
of inappropriate values of Z. These errors could be esti- 
mated using the approximate expression for the correction 
to the sound speed, namely, 



1 



dZ ), 



1 



3 + 5X - Z 



(7) 



T,p,X 

This error is found to be very small as compared to other 
uncertainties. It may be noted that if we take the par- 
tial derivative at constant Y instead of constant X, the 
right hand side of equation (Q) will be 6 times larger and 
the corresponding error will also be larger. In this work 
we therefore attempt to determine the X profile rather 
than Y profile from secondary inversion. With this choice 
the errors due to uncertainties in Z affecting the equation 
of state are negligible. Of course, the uncertainties in the 
equation of state itself will also affect the results, which 
can be estimated by repeating the inversion procedure 
with a different equation of state. Thus, while most of the 
results were obtained using the OPAL equation of state 
(Rogers et al. 1996), we have also done some inversions 



using the MHD equation of state (Dappen et al. 1988b 



Hummer & Mihalas [l988| Mihalas et al. |l988[ ). It turns out 
that the difference in these results is not significant and 
the computed luminosity decreases by 0.002L Q when the 
MHD equation of state is used instead of OPAL. It would 
seem that the secondary inversion results are not particu- 
larly sensitive to reasonable uncertainties in the equation 
of state. 

In principle, for any given Z profile it should be pos- 
sible to find T and X profiles which do not require any 
opacity modification as has been demonstrated by ST96. 
In fact, a choice of the parameter a = in equation (^J) 
produces a profile which requires very little opacity vari- 
ation and is similar to what is obtained if we were to use 
the procedure adopted by ST96. However, in this case it is 
not possible to ensure that the computed luminosity will 
necessarily match the observed solar luminosity Lq. This 
follows directly from the equations of thermal equilibrium 
(Eqs. (1,2)). Once the sound speed, density and Z profiles 
are known it is possible to integrate these equations to 
calculate the T and L r profiles, which will depend only on 
the central temperature T c . Clearly, by adjusting T c it is 
not possible to get both the correct luminosity and Y at 
the base of the convection zone. In fact, the solution is so 
sensitive to the choice of T c that a change in T c by mere 
1000 K, results in the value of Y in the convection zone 
to increase from 0.25 to a value greater than 1! Thus, it 
is not possible to alter the luminosity even by O.OOIL© by 
adjusting T c in the allowed range. Hence, the difference 
between the computed and the observed solar luminosity 
will give an estimate of uncertainties in primary inversions 



or the input Z profile, or the basic microphysics, such as 
the equation of state, opacities and the nuclear reaction 
rates. It may be difficult to separate out the contributions 
from each of these sources. We will try to examine this 
question in some detail in the following section. 

Of course, the inferred T and X profiles will depend 
on the assumed Z profile. In order to estimate this ef- 
fect, we try a few, different input Z profiles and the re- 
sults are shown in Fig. |^. It is clear that uncertainties 
in these profiles arising from those in Z are much larger 
than those due to other effects considered earlier. Further, 
these uncertainties decrease with r, because the effect of Z 
on opacities decreases as temperature increases because of 
increasing degree of ionization of heavy elements. There 
is a large uncertainty in regions immediately below the 
convection zone; however, in this region the value of Z is 
more reliably known from the measured value in the con- 
vection zone. Interestingly, the inverted temperature pro- 
file obtained using the Z profile PROF with Z sur f = 0.018 
comes close to that in Model S of Christensen-Dalsgaard 
et al. (1996). This is very likely because, with similar Z 
profiles, the opacities and hence the temperature gradient 
should be similar between two models. The variation in 
the sound speed is mainly reflected in the difference in X 
profile. The discordance between various profiles in Fig. 
could give an estimate of errors expected from reasonable 
errors in Z profile. 

The absolute X profiles as inferred using different Z 
profiles are shown in Fig. |6[ which also displays the profiles 
in some standard solar models with different treatment of 
diffusion. It is evident from this figure that the X profile 
just below the convection zone is much smoother than that 
in a standard solar model with conventi onal t reatment of 
diffusion (Christensen-Dalsgaard et al. 1996 ) suggesting 



that some turbulent mixin g probably takes place in this 
region (Richard et al . 1996 ). The X profile in Model 5 of 
Richard et al. ( 1996|) which includes turbulent diffusion 
is closer to the inverted profiles, though it appears to be 
shifted below the inverted profiles using Z sur f — 0.018, 
probably because it has larger Z sur { = 0.019 (and corre- 
spondingly higher Z/X). The shape of the AT-profile near 
the base of the convection zone is essentially independent 
of Z sm f, but depends on the actual profile used. Thus, if 
a flat Z profile like HOM or RICH is used the resulting X 
profile is also flat until about r = 0.68-R© indicating that 
this region is essentially mixed. However, if a Z profile 
with steep gradient near the base of the convection zone 
is used, then the resulting X profile also shows some weak 
gradient in that region. But in order to get a gradient as 
steep as that in the X profile of Model S, one requires 
a Z gradient which is about 5 times that in the profile 
PROF. Hence, if the X and Z profiles from similar treat- 
ment of diffusion in a solar model are used, the resulting 
profiles will not be consistent with helioseismic data un- 
less some process like turbulent diffusion is employed to 
reduce the gradients to zero at the base of the convection 
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Fig. 5. Relative difference in temperature and absolute difference in the hydrogen abundance X , between the Sun and 
the Model S of Christensen-Dalsgaard et al. ( 1996| ) as inferred using different Z profiles. The continuous line represents 
the results using the profile RICH, with Z smf = 0.018. The dotted (Z surf = 0.015), short-dashed (Z smi = 0.018) and 
long-dashed (Z SUI { — 0.02) lines shows the results obtained using the profile PROF, while the dot-dashed line displays 
the results using a homogeneous Z profile, with Z — 0.018. 



zone (Basu & Antia 1995 ; Basu 1997 ). We prefer to use 
inverted profiles with zero gradient in X or Z at the base 
of the convection zone for better accordance. 

Notice, around r — 0.25Rq the X profile in the Sun 
is steeper than that in the solar model. This can be seen 
more clearly from Fig. || which shows the difference in 
X profile between the Sun and a solar model. The steep 
positive gradient around r = O.25i? indicates that the X 
profile in the Sun is steeper than that in the model. 

In order to verify the seismically inferred composition 
profile, we have constructed a static solar model (Model 
INV) using the inverted profile for X as shown by the con- 
tinuous line in Fig. ^ and the model so computed is com- 
pared with results from primary inversions. This model 



also uses the OPAL equation of st ate (R ogers et al. 199G ) 
and opacities (Iglesias & Rogers 1996 ) and nuclear re- 
action rates from BP95, except for the pp reaction for 
which the cross-section estimated in the Section 4 is used. 
Since the surface hydrogen abundance X snl { and the mix- 
ing length parameter in these models are adjusted to get 
the correct radius and luminosity, the X profile has to be 
scaled by multiplying it by a constant factor and as such 
the resulting model does not have the same abundance 
profile as that inferred from inversion. Fig. [7] shows the rel- 
ative difference in sound speed and density between these 
models and the Sun, while the properties of this model 
are also summarized in Table ^. It is clear that the hump 
below the convection zone has more or less vanished. The 



Antia & Chitre: Seismic inversions for temperature and chemical composition 



9 






' 1 ' 






/ 












"**"*^ ^ / - 






/ 






















, , 1 1 , , 1 1 , 







0.4 



0.6 



r/R„ 



Fig. 6. The hydrogen abundance prohle as inferred using 
different Z profiles. The various line styles have the same 
representation as those in Fig. S. In addition the X prof iles 
in the Model S of Christensen-Dalsgaard et al. ( 1996| ) is 
shown by heavy dot long-dashed line, while that in the 
model 5 of Richard et al. (1996) is shown by the heavy 
short-dashed long-dashed line. 
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Fig. 7. Relative difference in sound speed and density be- 
tween the Sun and Model INV. The model INV has been 
constructed with the inverted profile of hydrogen abun- 
dance shown by the continuous line in Fig. ^. 



discrepancy in sound speed within the convection zone is 
likely to be due to uncertainties in the equation of state or 
the error in estimated radius, or may arise from er rors i n 
inversion due to influence of surface layers (Antia 1995|) . 
It can be seen that for this model the sound speed and 
density are very close to those in the Sun, more or less 
within the estimated errors in primary inversions. Some 
of the remaining differences could be attributed to errors 
in opacities, equation of state or nuclear reaction rates, 
which have not been adjusted while constructing these 
models or on account of errors in primary inversions. It 
may be noted that the neutrino fluxes in this model as 
well as other seismic models listed in Table 1 are signif- 
icantly lower than those in the standard solar model of 
BP95 with diffusion of helium and heavy elements, and 
marginally higher than those in the solar model of Turck- 
Chieze & Lopes (1993), which does not incorporate any 
diffusion of elements. 

4. Helioseismic estimate for the pp reaction rate 

The nuclear energy generation in the solar interior is 
mainly controlled by the pp reaction rate. The theo- 
retically estimated cross-section for this reaction varies 
from 3.89 x 1(T 25 MeV barns ( BP95) to 4.21 x 1CT 25 
MeV barns (Turck-Chieze & Lopes |1993| ). For convenience 
we denote the usually accepted value (Bahcall 1989| ) by 
So = 4.07 x 10~ 25 MeV barns, and express the cross- 



section in terms of So- Recently, there has been a claim 
that the pp nuclear reaction rate should be revised up- 
wards by a factor of 2.9 (Ivanov et al. 1997 ). Althoug h this 
claim has been contested (Bahcall & Kamionkowski 1997 ) 
on the nuclear physics grounds, it would be nice to have an 
indep enden t check from helioseismic data (Degl'Innocenti 
et al. 1997 ). The inverted profiles for T and X can be used 
to compute the total luminosity generated in the seismic 
models provided the nuclear reaction rates are assumed to 
be known. With the requirement that the Sun is in ther- 
mal equilibrium, the computed luminosity should agree 
with the observed value and that would impose some con- 
straint on the nuclear reaction rates. However, as we have 
seen in the previous section, the inverted profiles depend 
on the assumed profile for heavy element abundance. We 
therefore, investigate the effect of an assumed Z profile 
on integrated luminosity to constrain the nuclear reaction 
rate. 

Apart from the Z profile, there could also be uncer- 
tainties in the theoretically calculated values of opacities. 
In order to obtain constraints which are independent of 
errors in opacity we can consider X profiles, with coeffi- 
cients in equation (|^) chosen arbitrarily. These arbitrary 
profiles may not satisfy the equations of thermal equilib- 
rium with any reasonable estimate for opacities. However, 
using the inverted sound speed and assuming the X and Z 
profiles, it is possible to calculate the temperature profile 



10 



Antia & Chitre: Seismic inversions for temperature and chemical composition 



inside the Sun. Once these thermal and composition pro- 
files are known, the luminosity in corresponding seismic 
models can be computed. In order to estimate an upper 
limit on the pp nuclear reaction rate we can try to con- 
struct a profile which generates the minimum energy for 
the given sound speed and density profiles. Since the sound 
speed essentially constrains the value of T/n, where fi is 
the mean molecular weight of the solar material, it seems 
in order to cut down the energy generation one should re- 
duce T as well as /i to keep the ratio constant. It is clear 
that the minimum value of [i is achieved when X = 1 and 
Z = 0, i.e. when there is no helium or heavy elements 
present in the central region. From more detailed calcula- 
tion of energy generation rate, we have verified that this 
is indeed true, although strictly speaking, since the tem- 
perature is not high enough for helium burning reactions, 
the minimum energy generation occurs when X = 0, when 
there is no fuel to burn! But even a value of X = 0.005 
gives much higher energy generation rate as compared to 
X = 1 — Z in the core, because the temperature has to 
be increased when X decreases to keep the sound speed 
constant. Further, if the temperature is required to de- 
crease monotonically with radial distance, then even such 
profiles can be ruled out. Thus, leaving aside this unlikely 
possibility, the minimum energy generation occurs when 
X = 1 and Z = 0, when there is no helium in the core 
(Y = 0). 

For the case of a profile with X = 1 and Z = we can 
easily demonstrate that the computed luminosity in the 
resulting seismic model is about 0.617L Q when the usual 
nuclear reaction rates are adopted. Now if we increase the 
pp nuclear reaction rate for obtaining the correct solar 
luminosity with this profile, it turns out that the cross- 
section needs to be increased to about 1.625*0. It is clear 
that if the cross-section is increased beyond this value it 
is not possible to find any X profile (apart from the one 
where hydrogen is almost totally exhausted throughout 
the solar core), which will simultaneously yield the cor- 
rect sound speed and luminosity in seismic models. The 
exact limiting value of the cross-section will depend on the 
inverted sound speed and density profiles, but as we have 
seen in the previous section these uncertainties are very 
small. We can therefore, conclude that any value higher 
than 1. 655*0 is inadmissible even if arbitrary errors in opac- 
ities are allowed and the Sun is assumed to generate the 
observed luminosity. An increase in the pp nuclear reac- 



tion rate by a factor of 2.9 (Ivanov et al. 1997) is certainly 
ruled out by the helioseismic data. In fact, in actual prac- 
tice even the profile with Y = considered in obtaining 
this limit is unacceptable since one would expect signifi- 
cant amount of helium to be present in the solar core. If we 
consider a profile with Y = 0.2, which is still lower than 
the expected helium abundance, the limiting cross-section 
for the pp reaction drops to 1.27So. It is therefore evi- 
dent that, any significant increase in the pp cross-section 
is demonstably inconsistent with helioseismic constraints. 



We would like to add that there is no straightforward way 
to set a lower bound on this cross-section from such an 
analysis, as by increasing the helium abundance suitably 
it is possible to reproduce the solar luminosity even when 
this cross-section is significantly reduced, although such 
profiles may require inadmissibly large opacity modifica- 
tions. 

In the foregoing discussion we have allowed for arbi- 
trary errors in standard opacity tables. Even though such 
an analysis helps in illustrating that the helioseismic data 
are able to put severe constraints on nuclear reaction rates, 
the resulting bounds on cross-section are highly conserva- 
tive and are unlikely to be achieved in realistic situations. 
It would be possible to obtain more meaningful bounds if 
one allows only reasonable errors in opacities. There are 
two problems with this approach; first, it is difficult to de- 
fine what is a reasonable error in opacity and second, the 
error in opacities may have arbitrary variation with tem- 
perature and density, thus making it difficult to consider 
all possible variations even within the assumed limits. One 
possibility is to use the procedure outlined in the Section 2 
with a suitably large value of a in equation @ , to obtain 
the X profile which generates the correct luminosity and 
requires some minimum opacity variation for any specified 
nuclear reaction rates. We consider this approach later, 
but before that we adopt a simpler procedure by taking 
different Z profiles with a large range of Z sur f to see how 
the computed luminosity varies with Z. In this process, 
the opacity changes are accounted through changes in Z 
profiles. 



Using the Z profile with diffusion (Proffitt 1994) scaled 
to different values of Z sux { we can calculate the X profiles 
following the procedure outlined in Section 2, with a = 
in equation (|g). The total luminosity and neutrino fluxes 
in the resulting seismic models are shown in Fig. ||. It is 
clear that the integrated luminosity goes up with Z su[ { 
as a result of increase in opacities, but not very signifi- 
cantly - a variation of Z SUI { from to 0.06, results in an 
increase in the luminosity from O.87L to 1.13L Q . The 
range of Z values covered by these models is in all proba- 
bility much more than the expected uncertainties in the Z 
profile. With the allowance of a factor of two variation in 
Z sur { , one gets an error of about 5% in computed luminos- 
ity. An uncertainty by a factor of two is probably the most 
that is expected in Z , and hence we have only considered 
profiles where Z is scaled uniformly by the same factor. 
A slightly smaller change in Z in selected regions may 
also give rise to similar change in the resultant luminos- 
ity. From the results presented later where we try to adjust 
the Z profile to match the luminosity, it turns out that the 
required maximum change in Z is not much smaller than 
what is indicated by this simple analysis. For the purpose 
of this work we therefore estimate a reasonable error of 
5% in the luminosity arising from possible uncertainties 
in the heavy element abundance and/or opacities. Since 
this is much larger than the estimated uncertainties from 
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Fig. 8. The integrated luminosity and neutrino fluxes for 
seismic models as a function of Z sur {. All these quanti- 
ties are scaled with respect to the values in the standard 
solar model of BP95. The continuous line shows the com- 
puted luminosity, dotted line shows the neutrino flux in 
the chlorine experiment, short-dashed line shows the neu- 
trino flux in the Gallium experiment while the long-dashed 
line shows the 8 B neutrino flux. 

other sources we assume a total uncertainty of 5% in com- 
puted luminosity. There are uncertainties in other nuclear 
reaction rates which will also affect the computed lumi- 
nosity, but again if these are within limits given by BP95, 
the error in luminosity from these is only about 1-2%. 
Thus the integrated luminosity is consistent with the ob- 
served value within these uncertainties for a reasonable 
Z profile. It may be noted that all these results are ob- 
tained using the pp reaction cross-section to be Sq. If the 
recent value adopted by BP95 (0.9558So) was used, the 
computed luminosity would be about 4% lower, while for 
the normal value of Z the computed luminosity would be 
significantly lower than the observed value. This leads us 
to surmise that the cross-section for the pp-reaction rate 
needs to be increased to its earlier value given by Bahcall 
( |1989p . Similar conclusions were also reached earlier by 
Antia & Chitre (|i995|). 



In order to obtain a better estimate for the cross- 
section of pp reaction, we try to compute the luminosity 
using different values for the cross-section of the pp reac- 
tion, with the normal value of Z SUY f = 0.018. From these 
results we can identify the range of cross-section values 
which yield the computed luminosity within 5% of the 
observed value. This can be treated as the helioseismic 
estimate for the cross-section of pp reaction, which turns 
out to be (4.15 ± 0.25) x 10~ 25 MeV barns, where the 
quoted errors correspond to an uncertainty of 5% in the 



0.06 


1 


1 


I , 


0.04 








0.02 


■ - ' / 






o 













1 1 1 1 1 1 1 1 1 








0.02 


J' " N " 






0.04 








0.06 









0.2 



0.4 



0.6 



r/R e 



Fig. 9. The relative difference in temperature and abso- 
lute difference in sound speed betwe en the Sun and Model 
S of Christensen-Dalsgaard et al. ( 1996 ). These results 
have been obtained with the nuclear reaction rates from 
BP95 with the Z profile also adjusted to yield the correct 
luminosity. The continuous line shows the results obtained 
with pp reaction rate from BP95, while the dashed line 
shows the results obtained when pp reaction rate from 
Bahcall (1989) is used. The dotted lines denote the la 
error limits due to those in the input frequencies from 
GONG months 4-10 data. 



computed luminosity. This range is consistent with the 
value adopted by Bahcall ( |1989| ), but slightly larger than 
the more recent value adopted by BP95. 

All the inversion results presented so far were obtained 
using a = in equation (^]), which yield profiles that re- 
quire no opacity modifications, but the computed lumi- 
nosity may not match the observed value. It is possible to 
adjust the opacity or equivalently the heavy element abun- 
dance to obtain the correct observed luminosity by choos- 
ing a suitably large value of a in equation (^). However, 
such profiles may not be unique as only one parameter 
namely, the luminosity is fitted by adjusting the Z-profile 
in radiative interior. Nevertheless, we can obtain a possible 
solution which yields a seismic model with correct lumi- 
nosity. We use the nuclear reaction rates as adopted by 
BP95 (including that of pp reaction) for this study. Then 
the computed luminosity turns out to be about 0.94£q 
for the case of a — (i.e., no opacity modifications), but 
if we choose a large value of a, say 2500, the computed lu- 
minosity comes out to be close to the observed value. The 
resulting value of R K at each point can then be converted 
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Fig. 10. The Z-profile needed to obtain correct computed 
luminosity using the nuclear reaction rates adopted from 
BP95. The continuous line shows the required Z profile 
when pp reaction rate from BP95 is used, while the short 
dashed line shows the results obtained when pp reaction 
rate form Bahcall (1989) is used. The dotted lines denote 
the la error limits due to those in the input frequen- 
cies from GONG months 4-10. The long dashed and the 
dot dashed lines respectively, show the profiles PROF and 
RICH scaled to Z suli = 0.018. 



to equivalent variation in Z 7 and this gives the Z profile 
required to obtain the correct total solar luminosity. For 
this purpose we use the technique of simulated annealing 
to minimize the function F defined by equation (|J). Once 
the iteration is close to convergence we linearize about 
that solution and determine the actual X profile. The re- 
sulting T and X profiles are shown in Fig. [|, while Fig. [To] 
shows the required Z profile. These figures also show the 
estimated errors arising from those in frequencies. If the 
cross-section for the pp nuclear reaction is taken to be 
So then the required increase in integrated luminosity is 
much smaller and the resulting profiles are also shown in 
Figs. H and [To[ It is clear from these figures that with the 
BP95 reaction rates the profiles need to be modified signif- 
icantly to generate the extra luminosity and the resulting 
Z profile in the central regions does not look particularly 
realistic, as there is an unrealistic hump in the core. Thus, 
it appears that Z needs to be increased by about a factor 
of two to obtain the correct luminosity when the nuclear 
reaction rates as adopted by BP95 are used. If the Z pro- 
file is uniformly scaled by a constant factor (cf., Fig. ||) 
of two, there will be a similar change in the luminosity. 
Further, the estimated errors are also larger than those 
in Fig. [|, because the profiles are more sensitive to er- 



rors in primary inversions when the luminosity constraint 
is applied. When the pp nuclear reaction rate is adopted 
from Bahcall (1989) the resulting Z profile still has the 
same shape as before, but the height of the hump is much 
less and it is within 2a of the usual Z profiles including 
diffusion. 

The sensitivity of inverted profiles to nuclear reaction 
rates is probably due to helioseismic constraints which are 
being applied in this work. Since the density and sound 
speed are known from primary inversion, in order to main- 
tain the observed solar luminosity with the reduced nu- 
clear reaction rates, the temperature will need to be in- 
creased. The sound speed constraint from the primary in- 
versions fixes the ratio T / /x, where /x is the mean molecular 
weight. The mean molecular weight /x will also have to be 
increased, implying a decrease in the hydrogen abundance 
X. This will work against enhancing the energy produc- 
tion, and as a result, temperature will have to be increased 
substantially to keep up with the required nuclear energy 
production for maintaining the observed luminosity. If we 
allow for the departures from the inverted sound speed and 
density profiles then it may be possible to obtain the cor- 
rect luminosity without much modification in Z, but the 
resulting seismic model will not have the inferred sound 
speed in the core. 

5. Discussion and Conclusions 

In this work we have shown that the use of sound speed 
and density profiles obtained from primary inversions en- 
ables us to infer the temperature and hydrogen abundance 
profiles, provided the heavy element abundance profile as 
well as the microphysics like the equation of state, opac- 
ities and nuclear energy generation rates are known. The 
profiles so determined would define a seismic model, al- 
though the integrated luminosity in these seismic models 
may not necessarily match the observed value. The differ- 
ence may arise due to uncertainties in primary inversions, 
and/or the assumed Z profile, and/or the microphysics. 
While it is possible to adjust the Z profile to yield the 
observed luminosity, it is not clear if that is the correct 
solution, since the discrepancy may arise due to other rea- 
sons. We have attempted to estimate the extent to which 
various uncertainties can influence the luminosity to find 
that the effect of equation of state or primary inversions on 
computed luminosity to be fairly small. The dominant un- 
certainty arises from the nuclear reaction rates and opaci- 
ties (or equivalently the Z profile). It is difficult to separate 
out the influence of these two factors, but if we assume a 
reasonable error in one of these the other effect can be 
quantified. 

It turns out that if we use the nuclear reaction rates 
adopted by BP95, except for the pp reaction for which 
the older reaction rate from Bahcall ( 1989] ) is used, then 
the integrated luminosity with the normal Z profile is 
close to the observed value. It is thus tempting to con- 
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elude that these nuclear reaction rates, together with the 
current opacity tables and a Z profile including diffusion 
are consistent with helioseismic data. Similarly, from a de- 
tailed study of the base of the convection zone it appears 
that uncertainties in the current opacities at the base of 
the convection zone as well as the estimated Z/X values 
(Grevesse & Noels 1993) are fairly small (Basu & Antia 



1997). One expects opacities to be more reliably deter- 



mined in the solar core where temperatures are upwards 
of several million degrees. It is therefore reasonable to as- 
sume that there are no significant uncertainties in current 
OPAL opacity tables in the solar core. Of course, there 
could be some error on account of an inappropriate Z dif- 
fusion profile, but that is not expected to be too large. 

It is remarkable that, different values of cross-section 
for the pp-reaction have been adopted by various workers 
(Bahca ll |1989| ; BP95; Turck-Chieze & Lo pes |l993| ; Dar & 
Shaviv [1996 ) and recently Ivanov et al. ( 1997 ) have even 



suggested an increase in the pp reaction rate by a factor 
of 2.9. Since there is no experimental measurement of this 
cross-section, it would be interesting to indulge in an ex- 
ercise to estimate this cross-section helioseismically. From 
our results in the previous section it is clear that an in- 
crease in this reaction rate by a factor of 2.9 is essentially 
ruled out, even when arbitrary variations in opacities are 
allowed. The only X profiles which may yield the com- 
puted luminosity as low as the observed value with such 
nuclear reaction rates are those where hydrogen is almost 
completely exhausted (X < 0.002) in most of the core. 
In fact, even an increase by a factor of 1.65 in the pp 
reaction rate is inconsistent with helioseismic data, with 
no restriction on opacity. If the helium abundance is con- 
strained to a minimum of 0.2, then this limiting factor is 
decreased to 1.3. Thus we can firmly conclude that even a 
30% increase in the cross-section for pp reaction is incon- 
sistent with helioseismic data. However, these bounds are 
too conservative since unrestricted errors in opacity are 
permitted. 

Should we make the assumption, on the other hand 
that opacities are known to reasonable accuracy and that 
there is an uncertainty of up to a factor of two in Z SUI { , 
then this translates into an uncertainty of 5% in computed 
luminosity and the estimated value of the cross-section 
for the pp reaction turns out to be (4.15 ± 0.25) x 10~ 25 
MeV barns. This value is consistent with the estimate of 
4.07 x 10~ 25 MeV barns (Bahcall ^989] Dar & Shavi v [1996J ) 
or 4.21 x 10~ 25 MeV barns (Turck-Chieze & Lopes Jl993P , 
but slightly larger than the value of 3.89 x 10~ 25 MeV 
barns adopted by BP95. It thus appears that the estimate 
of the pp reaction cross-section adopted by BP95, needs 
to be increased by a few percent. With the adoption of 
the recent estimate of this cross-section, the Z profile will 
need to be modified by about a factor of two to obtain 
the correct computed luminosity. We cannot, of course, 
strictly rule out such Z profiles, but they appear unlikely 
to be realized in practice. 



The reliability of the inverted seismic profiles from the 
observed frequencies can be demonstrated by constructing 
a model (INV) with the inverted X profile. This model is 
found to be close to observations in many respects in- 
cluding the sound speed and density profiles through the 
solar interior. Thus, the base of the convection zone in 
model INV is found to be at a radial distance of O.7131i?0, 
which is consistent with the helioseismically estimated 



value (Christensen-Dalsgaard et al. 1991; Basu and Antia 



1997 ). Similarly, the helium abundance in the envelope 
of this model is found to be 0.2469, which is also close 
to th e helioseismically estimated value (Basu and Antia 
1995 ). This appears to indicate that the inverted X pro- 
file used in this model is close to that in the Sun. Note 
that this model has been constructed using the standard 
OPAL opacities and equation of state. It should be recog- 
nized that this model satisfying the seismic constraints is 
probably not unique, as it may be possible to construct dif- 
ferent solar models satisfying the helioseismic constraints 
by modifying the opacities or nuclear reaction rates or the 
Z profiles suitably. 

The inverted T and X profiles are found to be close 
to those in the Model S of Christensen-Dalsgaard et 



al. ( |1996[ ). It should be stressed that our technique is ab- 
solute in nature, and hence the fact that the resultant in- 
verted profiles are close to those of a standard solar model 
cannot be a coincidence, since the model profiles have not 
been used anywhere during the secondary inversion. Our 
results, therefore, appear to demonstrate that the temper- 
ature and hydrogen abundance profiles in the Sun are close 
to those in a standard solar model. The major noticeable 
difference arises just below the base of the convection zone, 
where the inverted X profile is smoother than that in the 
standard model. The X profile is in fact, sensibly flat in 
the region r > 0.68i?©. This is probably owing to some 
process involving turbulent diffusion just below the base 
of the convection zone, which is not accounted for in the 



usual treatment of diffusion (Richard et al. 1996 ). Such a 
mixing could smoothen the composition gradient and also 
explain the anomalously low lithium abundance in the so- 
lar photosphere. It should be stressed that the estimated 
uncertainties due to errors in the Z profile are fairly large 
and consequently, significance of the flatness of the profile 
may not be obvious. However, a mere increase or decrease 
in the opacity by a constant factor will not change the 
nature of the profile as similar results can be obtained for 
different values of Z sm f . Only if there is a sharp gradient in 
modified opacity over this narrow region (or equivalently 
a sharp gradient in the Z profile) it will be possible to 
obtain composition profiles which are not flat just below 
the convection zone. If the gradient in Z profile were to be 



increased by a factor of five over that in Proffitt (1994), it 



would be possible to get an X profile with gradient simi- 
lar to that in Model S at the base of the convection zone. 
Thus, composition profiles obtained using similar treat- 
ment of diffusion for both helium and heavy elements are 
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not consistent with inverted profiles unless the gradient 
vanishes as in the case of turbulent mixing (Richard et 
al. 1996). These results are consistent with conclusions 



drawn from the oscillatory signal in the frequencies (Basu 
& Antia 1994; Basu 1997), which also supports the pres- 



ence of turbulent mixing in this region. Similar evidence 
is also suggested by the inversion of sound speed (Gough 
et al. 1996[ ). All this seems to indicate that the region just 
bel ow th e convection zone is probably mixed (Richard et 
al. 1996 ) by some process. 

In contrast, in the central region around r = 0.25i?© 
the composition profile in the Sun appears to be steeper 
than that in the solar model, perhaps suggesting that mix- 
ing is unlikely to have occurred in this region of solar 
interior. From Fig. || it can be seen that 5X has a neg- 
ative gradient in the inner core around r — 0.1-Rq, which 
would imply that the X profile in the Sun is smoother 
than that in the model. This difference has presumably 
been considered as a hint of mixing in the core (Gough et 
al. 1996). However, considering the fact that the X gradi- 



ent is very steep in this region, the difference is extremely 
small and mixing if any, could only have taken place in the 
very early history of solar evolution or the mixing process 
is extremely slow. A more likely cause of this difference 
is the errors in nuclear reaction rates. It is also possible 
that this difference could arise from uncertainties in the 
primary inversion in the core. 

Using the inverted T and X profiles it is possible to es- 
timate the neutrino fluxes. From the results in Table 1, it 
appears that these neutrino fluxes are significantly lower 
than those in the standard solar model of BP95, with diffu- 
sion of helium and heavy elements. Some of the difference 
could be due to somewhat lower cross-section for pp reac- 
tion used by BP95. A part of the difference will also arise 
from the diffusion of heavy elements. As argued earlier 
there are good reasons to believe that the region imme- 
diately below the convection zone is mixed and hence the 
heavy element abundance will not increase as steeply as in 
the model of BP95. A reduction in Z value inside the core 
will reduce the opacities and hence the temperature and 
the corresponding neutrino fluxes. However, the computed 
neutrino fluxes in seismic models are significantly larger 
than the observ ed values. In fact, it has been found (Antia 
& Chitre 1997) that even if arbitrary variations in opac- 
ities are allowed it is not possible to reduce the neutrino 
fluxes in any two solar neutrino experiments simultane- 
ously to the observed values. Thus, it appears that the so- 
lution of solar neutrino problem should be sought in terms 
of neutrino properties, though the seismic models can be 
used to constrain these solutions. Since the neutrino fluxes 
in the standard solar model of BP95 are somewhat differ- 
ent from those in the seismic models, the constraints on 



the particle physics solution (e.g., Hata & Langacker 1997) 
could change when seismic models are used. 

We have demonstrated that our inversion technique 
produces reasonably well the thermal and composition 



profiles in the Sun's interior, with the knowledge of the 
sound speed and density inferred from the accurately ob- 
served frequencies, based on the mechanical and ther- 
mal equilibrium constraints governing the solar structure. 
These seismically determined temperature and hydrogen 
abundance profiles in the Sun turn out to be close to those 
obtained with a standard solar model. The small depar- 
tures could be due to a variety of processes arising from 
diffusion and uncertainties in nuclear reaction rates, equa- 
tion of state, heavy element abundances or even the pres- 
ence of a magnetic field. It is remarkable that the neutrino 
fluxes in the framework of the seismic model come out to 
be close to those predicted by the standard solar model, 
assuming that the opacities are not very different from the 
currently accepted OPAL values. There is thus a strong 
hint of the particle physics solution of the solar neutrino 
puzzle! 
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